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Abstrect 


An  addition  theorem  is  obtained  for  the  separated  solutions  of  the  reduced 
wave  equation  in  parabolic  cylinder  coordinates;  it  is  based  on  the  property  that 
the  reduced  wave  equation  is  invariant  under  translations  and  rotations  of  the 
coordinate  systeir..  An  interesting  feature  of  the  theorem  is  that  the  basic  set 
of  the  separsted  solutions  that  are  regular  and  one-valued  everywhere  is  split 
into  two  disjoint  subsets  such  that  the  addition  theorem  for  the  functions  of 
one  set  involves  functions  of  that  set  only. 

A  brief  discussion  is  devoted  to  the  problem  of  expanding  any  separated  one- 
vaJued  regular  solution  in  tenris  of  the  solutions  belonging  to  the  basic  set. 

Certain  integrals  which  arise  in  diffraction  problems  and  involve  functions 
of  the  parabolic  cylinder  are  evaluated  asymptotically  for  large  values  of  a 
parameter.  The  parameter  in  question  is  the  recipiocal  wavelength  of  the  inci- 
dent electromagnetic  field. 
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Introduction 

The  reduced  wave  equation     ^u  +  ku»OiE  invariant  under  translations 
and  rotations  of  the  coordinate   Fysteitio     WitMn  the  set  of  all  solutions,  the 
faiiiily  of  solutions  that  are  everywhere  regular  and  one-valued  must  form  an  in- 
variant subset  under  translations  and  rotations.     Now  suppose  that  in  a  certain 
coordinate   systeir     separation  of  the  variables  is  possible.     Then  a  translation 
or  rotation  applied  to  the  variables  of  a  separated  one-valued  and  regular  solu- 
tion leads  to  a   solution  which,   though  not  separated,  is  regular  and  one-valued 
again.     This  last  solution  must  therefore  be  expressible   as  a  lineai-  combination 
of  the  separated,  regular  and  one-valued  solutions  provided  that  these   are  com- 
plete in  the   space  of  the  regular  and  one-valued  solutions.     This  proves  the 
existence  of  addition  theorems  for  functions  arising  from  the   separation  of 
variables  for   Au+ku     «     0.     In  the  classical  case  of  polar  coordinates, 
the  resulting  regular  and  one-valued  functions  are,  in  two  dimensions,  the 
functions 

J^(lcp)  e^"^   ,  n  -  0,  ±  1,  t  2,    ...     . 

The  fact  that  the  non-regular  solutions  J  (kp)e  ^  and  HA  '^(kp)e  ^  also  have 
an  addition  theorem  can  be  derived  from  theorems  about  the  behavior  of  these 
functions  in  the  infinite  part  of  the  plane.  The  analogous  addition  theorems 
for  Mathieu  functions  and  spheroidal  wave  functions  have  been  derived  by  Meixner 
and  Schafke  ^  •' , 

In  the  case  of  the  coordinate  system  of  the  parabolic  cylinder,  the  sepa- 
rated solutions  that  are  everywl-iere  regular  and  one-valued  are  found  easily, 
but  their  completeness  shows  a  rather  curious  deficiency.  With  two  exceptions, 
all  plane  waves  can  be  expanded  in  convergent  series  of  these  regular  solutions. 
The  set  of  all  regular,  separated  solutions  splits  into  two  subsets  such  that 
the  expansion  of  a  plane  wave  involves  the  functions  of  one  set  only.  Consequent- 


ly,  the  addition  theorems  for  the  regular  solutions  are  such  that  they  involve 
functions  either  of  one  set  or  of  the  other  one. 

The  functions  of  the  parabolic  cylinder  have  previously  been  applied  tc  the 
diffraction  problem  arising  from  the  presence  of  an  electromagnetic  line  source 
in  the  interior  of  a  parabolic  cylinder  of  perfect  conductivity  (see  Magnus  L  -"and 
Buchholz  ^  ■• ) .  The  procedure  used  in  these  papers  has  been  criticized  by  Greenberg, 
et  al, '--'-»,  It  is  shown  in  the  present  paper  that  this  criticism  is  not  relevant. 
Also,  some  formulas  are  derived  which  describe  the  asymptotic  behavior  of  the 
diffracted  field  in  a  part  of  the  interior  of  the  cylinder  for  a  wavelength 
which  is  small  compared  with  the  focal  Jength  of  the  parabola. 

Part  I.  The  Addition  Theorems 

1,   Preliminaries 

We  shall  be  concerned  in  the  first  part  of  this  paper  with  the  problem  of 
addition  theorems  for  certain  functions  which  are  combinations  of  the  parabolic 
cylinder  functions. 

Generally  sf«aking,  we  consider  a  coordinate  system,  p,  C    for  which  the 
equation 

(1.1)  ^u  t-  k^u  «  0 

admits  separation  of  variables.  Then,  certain  solutions  of  (1.1)  will  be  of  the 
form 

(1.2)  u  -  co^(p,<r)  -  R^(p)S^(<r)  , 

where  v  is  a  complex  parameter.  In  general,  the  co  (p,cr)  are  regular  and  one- 
valued  only  for  a  subset  of  the  set  of  complex  numbers}  for  the  parabolic  cy- 
linder functions  this  subset  coincides  with  the  set  of  all  integers,  n. 

Suppose  now  that  the  x,y-plane  is  subjected  to  the  general  transformation 


(i.3) 


-  3  - 


x'     «       X  cos  a  +  y  sin  a  ♦  x       , 

o 


y'     «     -x  sin  a  +  y  cos  Cl  +  y       • 


Then     the  coordinates  p  -  p(x,y)   and  6"    ■   <5'(T.,y)  are  transformed  into 
p'   -  p(x«,y«)  and  <r«  -  C(x',y'). 
Write 

co^     -     M^(p   ,  <r  )  -  R^(p')  s^(«'). 

I 
One  of  our  principal  problems  will  be  to  express  these  co     in  terms  of  the 

w^CPjO"). 

Another  question  that  we  decide  is  that  of  the  conipleteness  of  the  set  of 
functions  oo     that  are  everywhere  one-valued  and  regular.     This  is  done  by  ob- 
serving that  all  but  two  plane  waves  can  be  expanded  in  terms  of  the  oo  .     It 
is  well  known  that  the  plane  waves  form  a  complete   system  in   the  sense  that  a 
regular,  one-valued  solution  of  (1.1)  can  be  approximated  in  any  finite  part 
of  the  plane  with  an  arbitrarily  Jsmall  error  term  by  a  linear  combination  of 
plane  waves. 

If  we  specialize   (1.3)  to  the  transformation 


(1.3)' 


X  ■       X    +    X 

0 


V. 


I 

y      =    y  +  y 


0 


we  obtain  a  formula  of  the  type 

(l.U)  CO  (p',6-')  .  J2    Y^"i  VP'^^  '^m^Po'^o^' 

yjhere  p     «  p(x  ,y  ),  (T   «  (TCx  ,y   ).     The  range  of  the   i, ,  m,n  is  the  group  of 

integers,   and  the  y.       are  numerical  factors. 

For  the  rotation 

»     X  cos  a  +  y  sin  a 
(1.3)" 


y'     ■  -X  sin  a  +  y  cos  a 


-  u  - 

we  obtain  a  relation  of  the  form 


(1.5)       co^(p',<r')  .  rA^"^(°)  '-JP*^')  » 

the  A    depending  only  on  a. 
m 


2.   Notations  and  auxiliary  formulas 
The  equation  which  concerns  us  is 

2 

(2.1)  ^u  +  ku-0,  u«  u(x,y), 

If  we  introduce  the  coordinates  of  the  parabolic  cylinder 

(2.2)  (K   +  i^)^  -  2(x  +  iy), 


or  eqvdvalently, 


X     .     |(C^  -  >|^), 
(2.3) 

y    "    ?  ■^ 

equation  (2.1)  becomes 

(2.U)  LE    +     L^    *     k2(42  ^     2j^     _     Q 

Putting 
(2.5)        u  -  f^(4)  fg^^^* 
we  obtain  for  f,  and  f_  the  following  equations: 

^  +  (w  +  kV)  f^  -  0, 


(2.6a)        '^^^ 


d4 

d  f  p  2  ? 

(2.6b)        =•  +  (-W  +  k'^yj^)  f,  -  0. 

^2  ^    '^ 
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Here  w  denotes  an  arbitrary  (complex)  constant. 

Particular  solutions  of  (2.6a)  and  (2,6b)   are 


(2.7a) 


f^U)  "  \  (^  (i*i)4)  ,  D^   (^yk   {l*t)K\  ,  D_^_^  U  VIT  (1+i)^    , 
C-v-l    f  ^  ^  ^^^^^^  ' 


(2.7b) 


^2^^  ^   "  °-v-l(>^  (l^i)^\     D_^^^  (-  yk  (i*i)v^y  D^(i  ylT  (l*i)v|), 
D^  (^-i/F  (l+i)^y 


where  v  »  - 


1      iw 
?"  5k  • 


The  functions  D  (z)  can  be  defined  by  a  power  series,  as  in    p.],  voi.   II, 


p.   117: 
(2.8) 


or 


(2.9) 


^  ,   ,       o(v-l)/2  -z^/U  -r  A      V    3       z     \ 


D  (z) 


2V/2,-^VU 


r(^) 


TSi  (      v       1        z     \       z 


r(-T) 


y?  r(-? 


1      V       3       z 


The  functions  -J  and  ^  are  explained  in  p.],  vol.  I,  chapter  6.  The  function 
^  is  sometimes  denoted  by 


■„  /     \   ,   a      a(a+l)  z 
^F^(ajc;z)  =  l  *  -•  z  +  ^-^^^  ^* 


.  •  •  • 


It  is  a  consequence  of  (2.9)  that  for  v-n-0,   1,   2,    ..., 


D  (z)     "     e 


-z  A  2-n/2   .  H^(2-V2,)^ 


-  6  - 
where  H  is  the  Hermite  polynomial  of  degree  n. 

Between  the  functions  D  and  D_  ,  there  exists  the  following  linear 
relation: 

(2.10)    D^(z)  -  '       e     C_v_i(iz)  +  e     ^-v-1^   M   * 

An  important  relation  which  is  a  consequence  of  (2.10)  is  obtained  by  putting 
V  ■  2n  and  replacing  z  l?y  iz: 

(-1)">^ 

^'•'''^       "7777"  "^Sn^i^^  =  ^.2n-i(-^^  ^  ^-2n-l(^^- 
(2n)l 

By  putting  v  »  2n+l  and  z  ■  -iy^  we  obtain 

(2.ab)   D3„^,(-ir>{)  ■  '"""^"*"'  |j-2n-2<Vl)  "  °.;n-2<-vV]  • 
These  two  relations  enable  us  to  write  the  functions 

(2.12)  V^'^)  =  \(Y4)B_y.l(rH)  +  ^v^"^^^^-v-l^"^^^* 

which  appear  in  subsequent  work,  in  a  simpler  form  when  v  is  an  integer.  Thus, 

^ir!^''^^^  '   ^2n^Y4)D_2„.-L(YV|)  *  T^2n^-riK2n.2.^"'V 

(2.13)  -  D2„(yO  [D.2r.-l(Yn)  *  I^-2n-i(-^^^] 

^^^^  D2„(r4)D,„(irn), 

since  D„  (z)  »  Do  (-z).  Similarly,  by  the  use  of  the  (2,11b)  we  find 
2n      tn 

(2.1U)        -2n.l(^>n)  -  ^-^X^^fT^    ^2n.l(^«^^2n.l(-n)- 
In  the  same  way  we  can  simplify  the  expressions  ^-^  ^  '• 


-  7  - 


and 


'-2n-l  -  ^-2n-l(Y4)D2„(Y>^)  *  ^.2n-l^-'^^^^2n^-r\) > 


It  will  be  necessary,  in  connection  with  the  expressions 


0>o_  »  CO, 


CO 


■2n*  "^2n+l»  "'-2n 


'^-2n-l*  ^°  consider  the  asymptotic  expansions  for  the  products  D  irk)!)       i(y^)> 
T^On^rO^jn^^^^^*  and  I52n+l^'''^^^2n+l^"^'^'^  ^*  ^°^  ■'"^^S®  "•  ^°^  ^^^^  purpose  we 


use 


(2.15)        D^(z)  ;ii  2"-^/^  exp  J  ^  log(-v)  -  j-  -  (-v) 


1/2 


where    jzl   is  bounded  and   |arg(-v)l  <  ^  ,   and   |v|  -s-ooj  we  use  also   (2.11a)   and 
(2.11b).     If  it  is  not  true  that    |arg(-v) |  <  j-  ,  then  we  must  first  use   (2.1C). 
From  these  relations  we  have  for    |v|   -j>oo: 


(2.16) 


D^(Y4)D_^,j^(YY|)as 


,-1 


VnTI 


exp^-Y(n+l)^/^V|+iC)  +  ^ 


+  exp^  Y(n+i)^/^(U->()  -  "^ 


V^^^V^^n^ 


(2n)l 


2  V^>/2mT 


exp^(2n+l)^/^Y(V^^^ 


+  exp^-(2n+l)^/^Y(i4+^) 


(2.17) 


+  expJ(2n+i)^/^Y(iC+>|) 


+  exp^(2n+])-'-/Vu-n) 


and 


-  8  - 


^2n.l^Y4)D2„,-L(-iYn) 


(2n-t-l)l 
2  >/2ir-/2H^ 


exp^-(2n+2)^/V(Y^+U)Uexp  J-(2n+2)^/^Y(U-^)!> 


(2.18) 


+exF  J(2n+2)^/^Y(i4-^  )  Uexp  J(2n+2)^/^Y(U+^  ) 


We  must  also  consider  the  asymptotic  expansion  of  the  product 


(2.19) 


\^rK)  t).v.i(Y>{) 


for 


V  =  -  X-  +  i^-,  X  real  and  tending  to  +oo.  For  X  -s>  oo. 


,-1 


nX       fX.in 


^.l/2.ix(Y«)^-l/2-ix(M>  ^  2"  exp^-  |  log  X  .  ^  -  ^r^rU  *V 


(2.20) 


i(^-X*   (^)^/^(4  *^))   I     • 


The  corresponding  expression  for  X  ->-oo  can  be  obtained  by  putting  X  »  -t 
and  letting  t  ->  +oo,  or  by  interchanging  ^  and  vi    and  letting  \.  ->oo  as  above. 


3,       Special  one-valued  solutions 

Before  entering  upon  a  discussion  of  the  addition  theorems  it  will  be  of 
some  value  to  remark  on  the  equation  of  one-valuedness. 

From  the  transformation  (2,3)  we  see  that  a  point  in  the   (x,y)-5ystem  has 
two  pairs  of  coordinates  in  the   (^,y^  ) -system,  namely  (C>>1  )  and  i-E,f-^)» 
Thus,  to  say  that  the  function  F(4, >y )  is  one-valued  in  the   (x,y)-plane  means 
that  F(^,>7)  »  F(-C,-v]).     Therefore,  a  product  fj^(YC)f2(Y">7)  is  a  one-valued 
function  in  the   (x,y)-plane  if  and  only  if  the   functions  f     and  fp  are  both 
even  or  both  odd. 

Given  f,(Y?)  and  f2(Y^)»  then  clearly 


-  9  - 

(3.1)  f^(Y4)f2(Y^)    +    f^i-r^)f^i'ryi) 

is  a  one-valued  solution  of  the  wave  equation,   although  one  in  which  the  variables 
are  not  completely  separated.     For  certain  purposes,  for  example  Cherry's  inversion 
formula  (  QL] ,  vol.   II,  p.   i2U),  these  are  the  proper  functions  to  use.     In  parti- 
cular, we  shall  need  the  functions 

(3.2)  co^(C,n)    -    \(ri^)^.^.i(r\)  *  ■D^irK)T>_^,^(-rV* 

They  are  fully  separable  for  integer  values  of  v,   according  to   (2.11)  and  (2.12). 

That  the  functions  (3.2)   form  a  basis  for  the  set  of  one-valued  regular 
solutions  of  the  wave  equation  follows  from  the  generating  function  relation  for 
the  parabolic  cylinder  functions 

^f         j.\  c+ioo  ^.  - 

F(x,y,t)  r-pr       I  .-v-1 

(3.3)  -74—    ^  sI^T^TV-'^''^"' 

where 

F(x,y,t)  .    -i4-  (-'  *  y')   -  ^  • 

U(i+t^)  i+t'^ 


Then,  if  in  (3.3)  we  replace  x  by  (T^  and  iy  by  Cr^     ,  where  0"  ■  -/k  (l+i),  we 

obtain 

c+ioo 


1+t        1+t 


^-ioo  yi+t 

2  G(4,)|,t), 

where  -i  <  c  <  0. 

If  we  introduce  an  angle  0  by  re  ans  of 

t  ■  tan  K-  , 

then  the  right-hand  side,  G,  of  (3.U)  is  a  plane  wave  incident  at  an  angle  0  to 
the  X-axis. 


-  10  - 
We  may  apply  the  Cauchy  residue  theorem  to  the  integral  in  (3.U),  which  has 
simple  poles  at  v  ■  0,  ±  1,  *  2,  . ..   .  If  we  close  the  path  of  integration  to 
the  left,  the  poles  inside  the  contour  are  v  *   -1,  -2,  ,,,,  and  we  obtain 

00 

(3.?)  G(4,>,  ,t)  =  V77^    y{-i}\{K,yi  H"""-^,        ( |t|  >  .!). 

n-1 

If  we  close  the  path  of  integration  to  the  right,  then  the  poles  inside  the  contour 
are  v  -0,  1,   2,    ,,.,   so  that 

00 

(3.6)  G(C,  v^  ,t)  «  vW?    Yl  (-l)Vn-l^"»  (1^1  <  1^- 

n=0 

In  (3.5)  and  (3.6)  the  inequalities  in  parentheses  indicate  the  range  for  vrhich 
the  series  converge. 

We  can  use  the  formulas  for  the  asymptotic  behavior  of  the  co  (<r^,C"v^)   for  large 
V  to  show  that  deformation  of  path  is  permitted.     Also,  we  can  prove  the  convergence 
of  (3.5)   and  (3.6)  by  means  of  (2.15)  and  (2.16).     Moreover,   (3.5)   and  (3.6)   show 
that  the  plane  waves  can  be  expanded,  •3xcept  for  those  which  are  defined  ty 
t  ■  i  1.     The  latter  are  the  waves  that  are  incident  at  an  angle  of  U5     to  the 
X-axis.     For  t  ■  +  1,  the  series  (3.5)  and  (3.6)  diverge. 


U,  Addition  theorem  for  translation 

From  the  considerations  of  the  previous  section  ve   are  led  to  the  addition 

theorem  for  translations. 

We  have 

'2    2     2   -2   ,  2, 


?'  ^'     ■  ^^    *  «o  ^o- 
From  the  definition  of  G(4,  vi  ,t)  we  see  that  the  following  relation  holds 
with  respect  to  the  above  translation: 


-  11  - 
(u.i)       G(c',n',t)    -  yCt^  •  G(4,v?,t) 

Now  let  us  put 


y[,t)   •  G(Cj,,  Y|  ^,t), 


v^ 


I    00 

4—   n    ' 


0 
00 


^~      n      ' 

0 


|t|  <  1 


Itl  >  1. 


Taking  coefficients  on  both  sides  of  (U.l)  gives  us 


00 

z 

o 


00 

z 

o 


(U.2) 


00 

z 

o 


.y^  •(^E(-i)"-.n.i(^o'^o)^") 


This  may  be  wrj  tten  as 


(U.3)    /^  (-i)"«.n-i^^  »^  )  -  ?  E  V"-c--i(^'"^Kt.i(€o»no)(--^) 

"n 

where  p,  6",  t',n  =  0,1,  2,  ...  . 

Hence,  the  addition  theorem  for  the  co   ,,  n  =  0,  1,  ...,  reads: 
'  -n-1 

p+  cf+t 


ff'  +  r 


The  coefficients  X    that  appear  above  are  defined  by 


X     »=  0       if  p  is  an  odd  integer; 


2r 


5-   \         {-lf{-  \)^ 
rl 


Equation  (U.U)  corresponds  to  the  case  |t|  <  1.  When  |t|  >  Ij 
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00 


vn       /v'    „  \x-n+l 


00 


^t-"*M.  v^T?  .  ((-1)"  03^(4, v^)f 


^'^•'^  .y^ff  (-1)"  0.^(4,, v^,)i 


oo 


-n+r 


Formula  (U.5)  provides  us  with  the  companion  addition  theorem  for  the  w  ,  n  "  0, 


f  I 


(U.6)  <^U  yV   )  'V^    JZ    \^'^^    "^cr^^ylKK^'^o^* 

"  p+a+t  ^ 


=n 


where  o,  cT,  t-O,  1,  2,    ...      . 

Formulas  (U.h)  and  (U.6)  exhibit  an  interesting  feature.     In   (U.U)  only  the 
functions  co      ,    are  related;  in   (U.5)  only  the  functions  co     are   related. 


5.       Addition  theoreni  for  rotation 

In  this  section  the  problem  is  to  obtain  an  addition  theorem  for  the  rotation 


(5.1)  i 


£       =     4  cos  S-    -    "^    sin  2- 


>l*  -    >^  cos  I    +4     sin  ^     . 


From  the  bilinear  continuous  generating  function  for  the  D^  (cf.    [1],  vol,   II, 
p.  126)  with  X  -  6E,,  iy  «  (fv^    ,  where  cr  »  -/k  (l+i),  we  get 


G(4,  yi  ,t)  '  —  exp  Jik     cos  0 

v/C7  ( 


2       2 

LZ^     +  sin  0  4>?1 


(5.2) 


"2^ 


c+ioo 


c-ioo 


-v-1 


4-7 r    CO    (<r   4,<rV}     )     dV 

sin(-viT)     v^     ^*      ^ 


The  trigonometric  functions  appearing  in  (5.2)  are  explained  by  the  substitution 
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t  »  tan  -  ,       frcm  which  i.t  follows  that 


(5.3)     < 


JL-t                    o,         2t 
jr     "  cos  (2,     jr 

1+t  1+t 


sin  |Z. 


Also 


(5.U)  G( 


Therefore 


I  2        2 

4*,  Y|*,t)  -  — i—  exp  Jik    Uos(9f-a)  ^-^4-  +  sin((?-a)^V|l 


/ — z^ 
(5.5)         G(«*,>)*,t)     .    121L:.    .  o(f,>,  ,t*). 


where 


.*       ^      /n      \         tan  0  -  tan  a 
t     «  tan(^-a)  " 

l+tan^l  tan  a 


If  we  put  tan  a  ■  A,   then  we  can  write  t     as 
(5.6) 


t*  .   (-A)    <   ^"^ 


rTAt  ;  • 


*  -K^ 


If  we  apply  the  residue  theorem  to  (5.2)  with  i>  V     replaced  by  4  ,    y? 
respectively,  we  obtain 


00 


G(c*,  v|*,t)  -  v^  r  i-if  co_^(c*,  N^*)t"-^ 


(5.7) 


AT? 


z:?^  "^ 


v-^  f<--.<M 


)t*-"-i . 


Putting  n-1  ■  m,  we  can  write 

(5.8)  X:  (-1)"'*^  CO  ^  ,(C*,  v^*)!"-  -  ^^-    Z  (-1)         "  ^  l(^->>?  )^       • 

yC7  ° 


Now 


we  must  expand  Vi+t*         /  Vl+t       and  t       on  the  right-hand  side   and 


-  lU  - 
then  compare  the  coefficients  of  similar  powers  of  t  on  both  sides.  We  shall 
then  have  the  required  addition  theorem. 
From  (5 •6)  we  have 

and,  for  m  >  0,  after  some  calculations: 

(5.10)    t   -  (-A)  }_  \^_{)\^i     2^1^"'^»~'"»^~""""'>"-^  )(-At)   . 


Here  ^F-   denotes  the  hypergeometric  series.      By  substituting   (5.9)   and  (5.10) 

in  (5.8)   and  comparing  coefficients  of  t"  on  both  sides,  we  obtain  the   following 

addition  theorer,  for  the  co       ,  : 

-m-i 

(^.11)  (-l)"*^  '^.m-1^^'''  "f^   '  "^    T:  (-1)'*^  '^-s-l(^'7  ^  ^6^^^* 

s»o 


wr.sre 


,m+s  ^       (s+n-Dl       ^  , „.-._   ....-2 


(5.12)         C  (A)  -  (-A)""*  ]Z      fSlTTnT    2V-"'-^i^-"-^'-*  ^• 

n=«o 

In  (5.12),  the  term 

(s+n-l)l 
(s-l)nl~ 

has  to  be  replaced  by  unity  if  s  «  n  «  0, 

It  can  be   shown  that  the  series  on  the  right-hand  side  of  (5.11)  converges 

for    |A|  <  1.     A  similar  addition  theorem  can  be  derived  for  the  (n^,  m  ■  0,   1,  2,  . 


6,       Expansion  problem  for  co 

We  novT  wish  to  investigate  the  possibility  of  obtaining  an  expansion  for 
the  one-valued  regular  solutions  oi     for  an  arbitrary  parameter  v  in  terms  of 
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the   solutions  (a,n"0,   ±1,   *2,    ...      , 
We  begin  id.th  the  formula 
c+ico 


^"^  •  2^  I    %(^*^>  ilET:^^^  '    G(4,v^,t) 


(6.i) 

•  exp  ^    ik 


Ll+t  1+t      -1 


(6.2) 


By  inversion  we  obtain  from  the  above  Mellin  integral 

<x^^{K,\)'iV^       G(?,Y^,t)  t^  dt 

0 

-  y^  (  I   G(f,vj  t)t^dt  +  J   G(C,>|  ,t)t^dt  j 

GO  N. 


1  ° 


where  -1  <  v  <  0.     Now,  if  an  expansion  of  the  kind  contemplated  were  possible, 
then  the  expression  in  parenthesis  in  the  last  member  of  equation  (6.2)  would 
be  equal  to 

0  o 

But  the  asymptotic  expansion  of  the  generic  tenm  of  each  series,  and  the 
application  of  the  Cauchy  ratio  test^ shows  that  the  expression  (6.3)  diverges. 
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Fart  II.     A  Diffraction  Problem 


1,       Formulation  of  the  problem 

We  consider  the  two-dimensional  wave  equation 


(1.1) 


9C  3>1 


in  the  coordinates  of  the  parabolic  cylinder.  We  are  interested  in  the  diffracted 
wave  u  due  to  a  source  at  a  point  on  the  axis  of  a  perfectly  conducting  parabolic 
cylinder.  Actually,  we  regard  the  source  as  a  line  source  generated  by  a  current 
in  a  wire  parallel  to  the  generators  of  the  cylinder. 
The  conditions  imposed  upon  u  are  the  following: 

(a)  On  the  parabolic  cylinder  the  values  of  u  are  the  negatives  of 
the  values  of  the  field  due  to  the  source;  i.e.. 


(1.2) 


(1.3) 


U(^,V|^) 


-H^^^kr), 


2^2 


^^■(M-^c^o)'*^-^. 


(1) . 


is  the  Hankel 


where  r"  »   (C"^   -£_>?_)"  +   ^"^  ~i  ' — ^—  ,   and  H^ 
function  of  the  first  kind. 

(b)  In  the  region  1^  |  <  |^_|>  the  function  u  should  satisfy  the 
wave  equation  and  should  be  regular,  "^  •  +  vi  represents  the 
boundary  of  the  parabolic  cylinder, 

(c)  The  function  u  should  satisfy  the  Soramerfeld  radiation  condition. 
The  correct  formulation  of  this  condition  for  our  problem  is  due 
to  Magnus'^-',     It   states,  in  the  case  k  ■  1,  that  for  .^  ->oo  and 
for  a  sufficiently  small  e  >  0, 


4-l|-iu< 


'>! 


,1-e 


H     dri 


'\ 
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where  u  is  the  complex  conjugate  of  u. 

A  function  u  satisfying  these  conditions  can  be  derived  from  a  formula  due 
to  Buchholz   (see    [2],  p.  177): 


H^^^kr) 


i  y?    1 


-d+ioo 


'^<5'-ioo 


n(-s.  J)    r(s.  ^)  yf  .   D.23.i/2(/=Hk  v^  ) 


(i.u) 


4sll/2(v'"^^)    •   °2s-l/2   (  V^^  V 


4s^l/2    (V=OT?o)    -^23-1/2   (V^^^) 


>  dB 


K    <  i 


where  D     are  the  parabolic  cylinder  functions  which  we  have  discussed  in  the 
first  part  of  the  paper,   and  the  E         are  defined  by 

„(0),„.  U/2     -2"/U     „     (V       1       z^\ 

From  (l.U)  we  can  write  for  the  function  u,  if  4  >  4, 

-(T+ioo 

"      V-ioo 


(1.5) 


•   ^28-1/2   (  ^''^  ^o^    •  ^-28-1/2   (  /^^   ^a> 


^-28-1/2   ^V^-^)       , 
.     -^ ds    . 


^:28.1/2(V^V 
In  formulas   (l.U)  and  (1.5),   s  =  <r+  if  .     In  (l.U),    |  (T  |  <  r  }  in   (1.5), 


0  <  6"  < 


c- 


It  remains  to  show  that  the  function  u  in  (1.5)  satisfies  the  conditions 
(a),  (b),  and  (c). 

For  the  boundary  condition  (a)  this  is  easLly  proved  by  putting  "'?  "  "^ 
in  (1,5).  For  the  investigation  of  the  other  conditions  it  will  be  more  con- 
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venient  to  confine  ourselves  to  the  case  in  which  E,     *    \     "0,  thereby  avoiding 
unnecessary  computational  difficulties.  The   expression  (l»5)  then  becomes 

<r  +ioo 

The  region  which  ve  consider  is  the  interior  of  the  parabola,  and  this  is  des- 
cribed by  the  inequalities 


-r|3  <  >?  <  v^ 


£  >  0. 


First  we  study  the  question  of  regularity  for  this  region.  This  involves 
the  properties  of  finiteness  and  differentiability  of  the  function  u«  The 
finiteness  is  a  consequence  of  the  behavior  of  the  expression 

(1.7)      D,,.,,,(«  y^  )  •  ^^ly'^'  •  r(-s  ♦  J)r(s .  ^) 

^-2s-l/2^^a^"'^'-^  ' 

for  large  s.  We  now  reqiire  asymptotic  expansions  of  the  functions  D  (z)  and 
E  '(z)  for  large  v.  From  [2],  page  kO,   formula  (2 7, a)  we  have 

Following  the  discussion  in  [l] ,  vol.  I,  page  279,  especially  formula  (12), 
we  have  for  |s|  ->cd: 

e21i/2(^^^  )-  ^eiM'/2^  p(l)e-i'^n'/2(.,)l/2(.,,^2^1/2 
(1.8)' 


•  y?7^  cos  r^  yiks  1 


We  choose  our  path  C  of  integration  in  the  right  half-plane  as  the  union  of  the 
two  paths: 
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2t  +  i  t^  0  <  T  <  00, 

-2r-ir^  -oo<t<0, 


For  f  <  0,  cos  v?  yiks   'v  cos  (T'+  i)  ^  "i^  .  The  absolute  value  of 
this  expression,  and  therefore  of  the  denominator,  is  bounded  away  from  zero, 
Thus,  for  T<  0,  the  quotient  of  the  two  E   -functions  is  bounded,  while 
for  f  >  0,  it  tends  to  zero  for  |  ^  |  <  ^  • 
As  for  D  (z),  we  have 

(1.9)  D„(^).2^/2(-l),-'/^^  (1-^,3,^). 


and 


From    QL],  vol.   I,  page  279,   (8),  with  v='26-^,a-^-^,     )t=s, 
X  "-ik^  }  we  obtain,   since  0  <  arg     >c  <  i  , 

°2s-l/2(^^''^  )  -^2^-3/^  .  e^^^^   .  4  >CTk     .  s-^/^^^  .   (-ik^^) 
.  e-=-(^/^)^^^'   .y?    cos    [ns  -   2 /Ili?  -  j]   . 
Now  s-^/^*^  =  s-^/^   .  e"  ^°^  ^     For    ^  >  0 

Re(s  log  s)  ~  -  J  r     +  U  t    log  r    . 
For    r  <  0, 

Re(E  log  s)  '^  -  5-  f  +  terms  of  order  Z"  or  f  log  |  "^  I  • 
Next  we  must  consider  the  factor 

cos  fns  -  2  yriiF    .  4  -  J.T  , 

c  V^  •  /^IT  '-'  i  r  c  v^    for     r  <  0. 

Hence  we  have,   apart  from  factors  whose  absolute  value  is  bounded; 

cos   Ins  -  2y^eikC     -  r  1 '^  e^P  •<   ^'^^  +  2|  r  |   C  7^  >  ,     for    t  <  Oj 


-   20  - 


and 


cos 


^s  -  2  /sikC^  -  jl 


Aj  exp  J    n   t       y         for      t   >     0. 


Now  consider    P(-s  +  p) 


PCs  +  J-);     we  have 


r(-s  +  ^)  PCs  +  J-)  -  2TI  exp  J   s  log  (-L.)  .  1  log  (-s^) 
If  In  s  >  0,   then  log  (^)   -  in;   if  Im  s  <  0,  then  log  (— )   =  -  in.     Therefore, 


r(-s  +  ^)  r(s+  ^)  --  2n  exp  Jis(sgn  Im  s)nl     •    (-s^)"^/^ 


(1.10) 


2ns-^/2  e-^^/'^       •  '^'^ 


(-  X"  +  2ti)n 


(-  t    -  2ir  )n 


if    r>  0 

if      t<  0 


In  absolute  value,  therefore. 


(1.11) 


r(-s  -r)r(s  -^  r)  I-  -^  •  e""'^  • 


Collecting  all  expressions  we  obtain  for   (1.7),  apart  from  factors  whose  absolute 
value  is  1: 

"2"       -rr^       "IT        Uriogltf        nt^     ^2r  .     -2t       -n/2Z^  ^'2s'l/2^^  ^^^^  ^ 
•   e  •s»a  ^'fi        •2        'e        -e'  i.  > '- 


(1.12)  2iTS       .  e    ^^  •   s^  •  e 


Toy 


^:2s.l/2(^a^^^^ 


which  tends  strongly  to  zero  as     Z    -ooo. 

Next,  we  study  the  property  of  differentiability  of  u  as  given  by  (1.6). 

p  o 

From  X  '  (C  -  ^    )/2  and  y  «  4  V|  ,   it  follows  that 


(1.13) 


1 
1 — ?  • 


9u 


3u      .    ^  8u 


i^r    -    Mir    +  ?  i^ 


9y 


Bf 


^^  /  4^  vj 
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There  is  no  problem  for  C  >  0.  But  for  4  =  0,^  really  means  tttt-  »  because, 

by  definition,  ^  >  0,     We  assume  rj    j'  0.  For  the  existence  and  continuity  of 

the  derivatives  5-  and  ^—  on  ^  «  0  it  is  necessary  that 
ax.  ay 

To  prove  this  we  observe  that 


3u 
34 


4-0     / 1        ^         ^       r(-s  *  h 


(1.15) 


ds. 


But  since  the  integrand  in  (l,l5)  has  no  poles  to  the  right  of  the  path  of 
integration,  condition  (l.lU)  is  proved. 

Lastly,  we  must  consider  the  radiation  condition.  Vie  shall  not  go  into 
detail.  We  merely  remark  that  we  can  base  our  work  on  a  formula  of  Magnus, 
(see  [5],  page  3U9,  formula  (19)).  This  formula  can  be  obtained  by  closing 
the  path  of  integration  with  respect  to  s  to  the  left.  We  then  obtain  an 
expansion  of 

H^^^(kr)  +  u 

in  terms  of  a  sum  of  products  of  the  type 

"28^-1/2  -^s^-1/^ 

where  the   s.     are  the   zeros  of  e\     ,  ,p ( Y|     '/-2ik   )  -  0  for  a  fixed     '^  p» 


2,       Remarks  on  the  case  of  short  wavelength 

In  this  section  we  are  interested  in  studying  the  behavior  of 


-  22  - 


(^•1)  e(°>    (V?  VI2IR  ) 

-2-^/2   ^^  f2L/2(^a>^^> 

for  large  k  «  2n/X..  This  requires  the  asymptotic  expansion  of  D  (z)  and 

F  (z)  for  large  z.  We  asaime  that  the  asymptotic  expansion  of  these  functions 

for  |z|  ->oo,  if  inserted  into  the  integral,  will  give  a  reasonable  approximation 

to  the  value  of  u  for  k  ->oo.  We  have,  for  |z|  -i>oo; 

2 
(2.2)  D^(z)  -    e""  /^   .  z^  , 


(2.3)  l^^"^^ 

2/,      2  -v+1/2    rc?) 

+  y7  •  e     '      •   (-J-) 


r(-^) 


Since  we  are  interested  in  the  case  in  which  the  line  source  lies  on  the  axis 
of  the  parabola,  we  take     vi  -  0 .     Then  we  have 


'^"'  e.[-  i|£je.ri»<=-  ^']f^  (W,-A«.p-|l!]  -|-^  ,...V-V- 


5 ;=  :nji^i 


.p[.  ^Je.[-M..  ^.]^  <..,^-/^e.ft^]^-.-.^r 
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We  close  the  path  of  integration  to  the  right  and  apply  the  residue  theorem 
(simple  poles  at  s  =  r-+  n,  n  =  0,  1,  2,  ...)»  obtaining 

(2.5)  u^       \         exp     ^  {I     -  I   )\    ; Y-       ■^     » 

"  L*^         °  J     y?  >/l2ik     n=n  Pfn+n 


V|  V^Tsik     n=0  r(n+l) 


•«o^f        ^o^ 


where  w  >»  (    — 5-  \    ;  — 5-  ■  coordinate  of  the  cusp  of  the  caustic  of  reflected 


a '  ^a 

rays .     Since 


n«=0 
the  expression  (2.$)  becomes 


(3.6,  .^Ul±      ''^- 


C  ^  -" .  ,.  -  «)V., 


^4  -«')" 


^. /:^iif 


(1  -  .)-^/2 


la 


On  the  axis  of  the  parabola,  between  the  cusp  of  the  caustic  (in  the  limiting 
case  of  geometrical  optics)  and  the  vertex  of  the  parabola^ the  field  is  everywhere 
of  the  same  order  of  magnitude  with  resyject  tc  k.  The  fact  that  our  formula 
breaks  down  for  w  =  1  indicates  that  this  order  of  magnitude  changes  at  the  cusp. 

In  the  integral  (2.1),  the  behavior  cf  D  (z)  and  E   (z)  should  be  studied 
for  large  parameter  v  as  well  as  for  large  |z|.  However,  the  results  of  the  pre- 
ceding section  suggest  that  the  contribution  for  large  v  may  be  negligible. 
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